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Abstract. We study the cosmological constraints on active-sterile neutrino 
oscillations v e <-> v a for the case when v s is partially filled initially, i.e. < 5N S < 
1. We provide numerical analysis of the cosmological production of 4 He, in the 
presence of v e <-> v a oscillations, effective after neutrino decoupling, accounting for all 
known oscillations effects on cosmological nucleosynthesis. Cosmological constraints 
on oscillation parameters corresponding to higher than 5% ^He overproduction and 
different non-zero initial populations of the sterile state 5N S < 1 are calculated. 
These generalized cosmological constraints corresponding to 5N S > are relaxed in 
comparison with the SN S = case and the relaxation is proportional to 8N S . 
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1. Introduction 

Atmospheric, solar and terrestrial neutrino oscillations experiments provided evidence 
for neutrino oscillations. The atmospheric neutrino data analysis confirmed the <-> v T 
channel as the dominant atmospheric oscillation solution. The results of the terrestrial 
experiment KamLAND and the results of SNO salt phase prefer flavor oscillation LMA 
solutions to the v e «-> u s one for the solar neutrino problem. Active-sterile neutrino 
oscillations are discussed as a supplementary sub-dominant channel. Even after the 
confirmation of LMA solution as a dominant solution of the solar neutrino problem, 
the oscillation parameters are not precisely known - the completeness of LMA solution 
was questioned and the scope for some possible sub-dominant transitions was explored 
[HElElSEliniQIE]- The available results of the neutrino oscillations from the solar 
experiments, including the later SNO ES, CC and NC data, slightly favor the existence 
of a small sterile oscillating sector. Thus, it is still interesting to explore the cosmological 
influence of sterile oscillating neutrinos. 

Neutrino v e <-> v s oscillations affect the expansion rate and the neutrino involved 
processes in the early Universe and in particular, influence Big Bang Nucleosynthesis 
(BBN) This allows to put stringent constraints on neutrino v e <-» v s oscillations 

from BBN considerations. The whole LMA and partially LOW v e <-+ v s solar oscillation 
solutions and atmospheric v e «-> v s solution were excluded by cosmological considerations 

nn na H3i nn na iiei nzi na nni 1201 1211 1221 12^1 1211 

Detail reviews discussing neutrino role in cosmology and neutrino masses and 
mixings may be found in refs. [213 123 EE EH] • 

Cosmological constraints were obtained usually assuming zero sterile neutrino state 
population before neutrino oscillations epoch, 5N S = p u Jpt 9 e — 0, where pf% is the 
equilibrium electron neutrino energy density. However, sterile neutrinos production 
is predicted by different types of models, like GUT models, models with large extra 
dimensions, many-fold Universe models, mirror matter models, neutrino oscillations 
models, etc., and may be present at the onset of BBN epoch. Hence, the degree 
of population of v s and its initial energy spectrum distribution depends on the v s 
production model, and in general SN S > 0. 

The general case of initially (before oscillations) non-zero sterile neutrino population 
was examined and its influence on the oscillation effects at BBN epoch was discussed 
in reference [21] • It was found that the kinetic effect of oscillations is sensitive to 
the initial population of v s . BBN constraints on neutrino oscillations parameters in 
the general case of sterile neutrino state, partially occupied initially, were discussed 
in reference [2S] for oscillations effective after electron neutrino decoupling. The case 
of 3% "He overproduction was considered. The presence of non-empty sterile neutrino 
< 5N S < 0.54 before oscillations was shown to lead to strengthening of the cosmological 
constraints proportional to the SN S value. 

Recent studies of the systematic uncertainties of primordial Tie abundance suggest 
as a more reliable value of Tie observational uncertainty 5Y p /Y p > 5%. In what 
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follows we derive 5.2% cosmological constraints, corresponding to SY p /Y p > 5% 'Tie 
overproduction on u e <-» v s oscillation parameters for the same type of non-equilibrium 
neutrino oscillations. We show that in contrast to 3% Tie constraints, in the 5Y P /Y P > 
5% case the non-zero initial sterile population leads to a relaxation of the constraints on 
oscillation parameters proportional to the value of SN S , and finally to their alleviation 
for 5N S = 1. 

In the next section, we discuss briefly the primordial Tie production in the case of 
non-zero sterile population at the start of oscillations. In the last section we present 
the cosmological constraints corresponding to Tie-overproduction 5Y p /Y p > 5% and 
different levels of initial population of the sterile state and compare them to SN S = 
case. 

2. Production of primordial "lie and neutrino oscillations 

2.1. Primordial 4 He 

"lie is the most abundantly produced, most precisely measured and calculated element 
among the primordially formed light elements. It has also a simple post-BBN chemical 
evolution - it is only produced in stars. Therefore, it is the preferred light element used 
for obtaining limits on nonstandard physics. Primordially produced Tie is calculated 
with great precision EDI EH E21 EH! = 0.2485 ± 0.0005. The predicted Tie value is 
in agreement with the observational data for Tie inferred from astrophysical observations 
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Currently there are significantly different observational determinations of the 
primordial 4 He abundance, which point to the existence of greater systematic errors 
than assumed before. Izotov and Thuan have found Y p « 0.242 ± 0.009 for 82 HII 
regions and Y p = 0.2421 ± 0.0021 for 7 preferred HII regions, and Y p = 0.2429 ± 0.0009 
for an extended data set of 89 HII regions |35] 136] . Olive and Skillman have re-analyzed 
their data and gave Y p « 0.245 ± 0.013 for all 82 HII regions and Y p = 0.2491 ± 0.0091 
for the 7 preferred targets (HZj- The real uncertainty on the primordial Tie abundance 
may be larger than previously assumed 3%. Present determinations of primordial Tie 
abundance from observations of extragalactic HII regions, indicate a significantly grater 
uncertainty for the Tie mass fraction EH dHl HHj namely dY/Y ~ 5%. | 

Thus a derivation of cosmological constraints on oscillation parameters correspond- 
ing to about 5% "He overproduction seems relevant. 

The primordial Tie abundance essentially depends on the freeze-out of nucleons, 
which occurs when in the process of expansion the weak processes rates T w , governing 
the neutron to proton transitions, become comparable to the expansion rate H(t). So, 
the primordially produced mass fraction of Tie Y p ~ 2(n/p)f/(l + n/p)f, is a strong 
function of relativistic degrees of freedom at BBN epoch, which enter through H(t), 

|| See also the pioneer papers discussing the possibility for higher systematic errors in the determination 
of primordial ^Re abundance, which was believed to be much smaller than at present |41l 1421 135] . 
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where H(t) ~ ^G^g e ff T 2 . § Y p also depends on the electron neutrino characteristics, 
namely neutrino energy spectrum, number densities N v and the neutrino-antineutrino 
asymmetry, entering through 

2.2. Oscillations effects 

Non-empty initial sterile neutrino population, present at the start of oscillations, 
influences BBN by (i) increasing the expansion rate and (ii) suppressing the kinetic 
effects of v e <-> v s on BBN. 

(i) The presence of 5N S means increased effective number of neutrino species in 

1 /2 

equilibrium during BBN, and leads to faster expansion of the Universe, H(t) ~ gj^ 
causing earlier n/p-freezing, Tf ~ (gef/) 1 / 6 , hence, ^Re is overproduced. This 
dynamical effect is parametrized through SY^ ~ 0.0135iV s . 

(ii) In case of oscillations between non-equilibrium sterile neutrino state < 5N S < 1 
and electron neutrino, proceeding after v decoupling, v e <-> v s oscillations lead to 
considerable and continuous deviations from the equilibrium v e spectrum (spectrum 
distortion) and production of neutrino-antineutrino asymmetry. || 

This effects nucleons kinetics in the pre-BBN epoch JTOl HH Elj This kinetic 
effect may be parametrized in terms of additional neutrino 5N kin and will be 
denoted further SYki n ~ 0.0135Nki n - Non-zero 5N S suppresses the kinetic effect 
of oscillations: SN^ decreases with the increase of SN S , hence 5N S leads to a 
decrease of the overproduction of^ie. 

For a wide range of oscillation parameters the kinetic effect of oscillations is large 
during the period of freezing of the nucleons and, therefore, affects BBN. This effect 
is strong even when there is a considerable population of the sterile neutrino state at 
the onset of the electron-sterile oscillations. The v e energy spectrum distortion is the 
greatest, if the sterile state is empty at the start of oscillations, SN S = and decreases 
with the increase of the sterile state degree of population [2JJ . The same behavior is to 
be expected for the kinetic effect of oscillations. 

The total effect of (i) and (ii) can be approximately described by SY p ~ 0.013<5iV 
SN = SN S + SN kin , where SN kin = 6N££ X (1 - SN S ) and 5N%™ is the kinetic oscillations 
effect, corresponding to 5N S = 0. The expression presents a good approximation to the 
numerically calculated dependence of the kinetic effect on the initial population of u s , 
derived in reference |27j : 

§ Due to its strong dependence on g e ff, 4 Hc abundance was used to put constraints on the number of the 
relativistic particles during BBN epoch |45l I4HI Wf\ . parametrized by 6N U = N u — 3. Cosmological data 
(except BBN) from CMB, LSS, SN, are not very restrictive on the extra light particles (4*%1 14*31 15T] . 
|| Neutrino-antineutrino asymmetry is generated during the v e <-> v s resonant oscillations |14| . This 
dynamically produced asymmetry suppresses oscillations at small mixing angles. Therefore, it leads 
to underproduction of 4 He and alleviates BBN constraints on oscillation parameters compared to the 
case without the account of asymmetry growth. We have accounted in this work for its effect as well, 
however the effect is sub-dominant and will not be discussed further on. 
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Figure 1. The solid curve presents the frozen neutron density relative to nucleons 
X n = n^/n nuc as a function of the sterile neutrino initial population, at 5m = 1CP 7 
eV 2 and sin 2 26 = 1. The dotted curve presents the kinetic effect, while the dashed 
curve presents energy density increase effect. 



There is an interesting interplay between the different effects which 5N S ^ exerts 
on oscillations and on BBN with oscillations and hence, the 5Y p production and the 
cosmological constraints on oscillations parameters for the case SN S ^ differ from the 
ones derived in references [TBI EH for SN S = 0. 

As found in [2EJ, for SY P /Y P > 5% corresponding to SN^ X > 1, the suppression 
effect (ii) dominates over the dynamical effect (i) of 5N S ^ 0. Hence, the total 
effect is a decreasing function of SN S , i.e. "He overproduction decreases with SN S 
(in comparison with the case SN S = 0) and correspondingly the BBN constraints on 
oscillation parameters relax. In the opposite case SY P /Y P < 5%, corresponding to 
SN^ X < 1, the dynamical effect (i) dominates and the total effect is increasing with 5N S . 
I.e. ^e overproduction increases and the BBN constraints on oscillations strengthen in 
comparison with the case 5N S = 0. In the case SY P /Y P = 5% the constraints for 8N S ^ 
coincide with the ones for SN S = 0, due to the cancellation of the two effects (i) and (ii). 

For illustration of SN^ X > 1 case, which we consider further on, we present in 
Figure 1 the effects (i) and (ii) on 4 He overproduction at 5m = 10~ 7 eV 2 and sin 2 29 = 1. 
We have studied numerically the contribution of these effects on neutrons to nucleons 
freezing ratio X n = n^/nnuc for different 5N S . The primordial yield of helium to a good 
approximation is expressed through X n : Y p ~ X n exp(— t/r n ), where r n is the neutron 
lifetime. For the chosen set of parameters the ^e overproduction decreases with the 
increase of 5N S . The suppression effect (ii) of 5N S dominates and a relaxation of the 
cosmological constraints compared to SN S = must be expected. 
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In the next section we present and discuss the calculated cosmological constraints 
on oscillation parameters corresponding to higher than 5% ^le overproduction and 
< 5N S < 1. 

3. Cosmological constraints on oscillation parameters 

BBN constraints corresponding to SY p /Y p = 3% overproduction of ^He and non-zero 
initial population of the sterile neutrino 5N S < 0.54 were calculated recently [2*S] . As 
far as 5Y p /Y p = 3% corresponds to SN^ X < 1, the constraints strengthen with the 
increase of the SN S value. They increase the BBN 3% ^e exclusion region for oscillation 
parameters corresponding to 5N S = towards smaller 8m 2 . 

Having in mind the existence of a large systematic error of "Tie measurements, in 
this work we provide a numerical calculation of cosmological constraints corresponding 
to 5Y P /Y P > 5% "He overproduction (8N^ X > 1) and different initial degrees of sterile 
neutrino population in the range < 5N S < 1. We have chosen SY p /Y p = 5.2%, i.e. 
a value slightly higher than the critical one 5%, inorder to illustrate the possibility for 
relaxation of the cosmological constraints on oscillations for SN S ^ 0. 

Our numerical analysis has shown that cosmological constraints corresponding to 
5.2% ^He overproduction relax with the increase of d~N s values. Up to 5N S = 0.5 the 
cosmological constraints are slightly relaxed in comparison with the case SN S = 0, 
however, for higher 5N S values, the constraints relax noticeably. The reason for that is 
the predominance of (ii) effect (the suppression of the oscillations kinetic effects) over 
the dynamical effect (i) for the given uncertainty of ^e. All cosmological constraints 
corresponding to SY p /Y p > 5% will have such behaviour, namely, they will be relaxed in 
comparison to the constraints for 5N S = 0. (Vice versa, the constraints corresponding to 
lower than 5% ^e uncertainty will be more stringent than the ones for initially empty 
sterile neutrino state.) 

In Figure 2 SY P /Y P = 5.2% BBN constraints are presented for different values of 
the initial population of the sterile state, namely the lowest dashed contour corresponds 
to a zero population, the solid curve corresponds to 6N S = 0.5, the dotted and the dot 
dashed contours to 5N S = 0.7 and to 6~N S = 0.9, respectively. 

Another interesting result is that there are considerable constraints even for a very 
high 5N S values for that really high ^e uncertainty. I.e. the constraints on neutrino 
mixing parameter vanish only when the sterile state is in equilibrium before oscillations, 
when the kinetic effect due to neutrino spectrum distortion disappears, i.e. SN S = 1. 

4. Conclusions 

We have studied BBN constraints on neutrino v e <-> v s oscillations for the specific 
case when the sterile neutrino is partially filled initially < 5N S < 1. Non-zero SN S 
has two-fold effect on BBN with neutrino oscillations: a dynamical effect leading to 
overproduction of He-4 and a kinetic effect, leading to underproduction of He-4 in 
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Figure 2. The dashed contour presents SY V /Y P = 5.2% BBN constraints for 6N S = 0, 
the solid curve corresponds to 5N S = 0.5, the dotted and the dot dashed contours — 
5N S = 0.7 and to SN S — 0.9, respectively. The resonant oscillations case corresponds 
to l.h.s of the figure, the non-resonant one to the r.h.s.. 



comparison with the case of deltaN s = 0. So, depending on the interplay between these 
opposite effects, the cosmological constraints may be either relaxed or strengthened. 

We have provided detail numerical analysis of the BBN production of ^e, Y p , in the 
presence of v e «-» v s neutrino oscillations, effective after electron neutrino decoupling. 
We have calculated and discussed cosmological constraints on oscillation parameters, 
corresponding to higher than 5% uncertainty of helium-4, for non-zero initial population 
of the sterile state 8N S < 1. It was found that the cosmological constraints on oscillation 
parameters relax with the increase of SN S . The cosmological constraints for the cases 
SN S < 0.5 are slightly changed in comparison to SN S = case, however, for bigger SN S 
the constraints are relaxed considerably and for 5N S = 1 they are alleviated. 

Resuming the results of the works discussing 5N S effect on BBN: Cosmological 
constraints corresponding to higher than 5% uncertainty of helium-4, relax with the 
increase of the initial population of the sterile state, while the constraints corresponding 
to lower than 5% uncertainty of helium-4 strengthen with 5N S . 

It is remarkable, that in case of BBN with non-equilibrium oscillations between 
electron and sterile neutrinos, it is possible to obtain cosmological constraints on 
oscillation parameters even in the case when the ^He abundance is known with 
uncertainty greater than 5% (Actually, it is possible to derive constraints on neutrino 
oscillation parameters for He-4 uncertainty up to 32% in the resonant oscillations case, 
and up to 14% uncertainty in the nonresonant oscillations far as these are 

the maximal possible helium overproduction values [12]). The cosmological constraints 
persist while initially the sterile state is non equilibrium. When SN S = 1 initially, 
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the kinetic effect of oscillations is zero and it is not possible to obtain constraints on 
oscillation parameters in the discussed BBN models with non-equilibrium electron-sterile 
neutrino oscillations, as well as in BBN models with equilibrium electron-sterile neutrino 
oscillations. 
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